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LINEARIZED  THEORY  OF  FOIL  IN  SUBSONIC  FLOW 


G.  A.  Dombro vsk iy  and  V.  S.  Tkalich 


Institute  of  Hydromechanics,  Academy  of  Sciences,  UkSSR,  Kiev 


The  main  results  of  the  linearized  theory  cf  a subsonic  gas  flow 
past  a foil  is  contained  in  the  foraula 


where  is  the  pressure  coefficient  at  points  on  the  foil  in  the  gas 
flow,  c«  - pressure  coefficient  at  corresponding  points  on  the  sane 
foil  in  a flow  of  an  incompressible  fluid,  - value  of  the  Bach 
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number  at  an  infinitely  distant  point  in  the  gas  flow. 

As  our  first  approximation  formula  (0.  1)  is  correct  for  the 
arbitrary  relationship  between  density  p and  pressure  p.  In  the 
present  study  we  demonstrate  that  for  each  Hoe><  1 we  can  find  a 
dependence  between  p and  p at  which  formula  (0.1)  becomes  precise. 
Because  of  this  result  it  is  possible  to  approach  evaluation  of  the 
accuracy  of  the  linearized  theory  differently  and  to  compose  a more 
graphic  linearized  theory  with  methods  based  on  approximations  of  the 
indicated  connection  between  density  and  pressure  and,  primarily,  by 
the  well  known  approximte  method  of  Chaplygin. 


§ 1 


Let  us  proceed  from  a linear  system  in  hodcgraph  variables, 
reduced  to  the  symmetrical  form  of 

<i.i) 

where  * is  the  velocity  potential,  <|>  - current  function,  9 - angle  of 
slope  of  velocity  vector,  s (v)  - function  of  velocity  modulus  v,  K (s) 
- Chaplygin  function  [1,  3]. 
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In  his  tine  B.  Sauer  [2]  examined  the  case  of 

<1.2)  K = n*o\  o = *-*„ 

where  n and  s0  are  actual  arbitrary  constants. 

There  was  a similar  study  of  this  case  in  monograph  [ 1], 
According  to  [1],  the  general  solutions  of  system  (1.1)  with  the 
Chaplygin  function  (1.2)  and  corresponding  formulas  for  velocity  v, 
density  p,  and  pressure  p are  written: 

(1-3)  ? = n*9le[^(a)-jF(j)dj]f  = 

where  F(b)  represents  the  arbitrary  analytical  function  of  complex 
variable  » = * — 


0-4) 


Bo 


P(°)  = 


A sli  o -f  ch  o * 

A sh  a + ch  o 


n*o  (.A  (sh  a — a ch  a)  + (ch  a — « Bh  «))  * 
_,_v  Bt  1 . „ 

~ nU  TthT+1  + C' 


where  A,  B,  C are  arbitrary  constants. 


Also  interesting  is  the  dependence  n = H («)  which  can  be 
determined  using  (1.2)  and  the  second  of  the  relationships  of 


(1-4) 


f ' 

r 
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according  to  formula  [3] 


(1.5) 


M - VI””*? 


The  transition  to  the  physical  plane  of  the  gas  flow  Z = X ♦ iY 
is  achieved  through  formula  ([1],  page  63) 


(1.6) 


where  V9  = e*m 


Significant  for  us  in  the  present  study  is  the  particular  case 
of  A = -1.  Assuming  that  A = -1,  we  get 


(1.7) 


v(o)  = Boe\  p(o)=^JTTT, 

„ *****  a,  /.V  __  v'-jT+T 

P(3)  — C\  Jut*  M(o)  — ■ a-jj_  J , 


where  Ct  * C - B*/2n*. 

It  is  significant  here  that  the  formula  of  transition  to  the 
physical  plane  is  also  simplified: 
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(1.8) 


z=-r{-£['’  <»>-f  <>H+4  J-J  ■£''«}• 


5 2 


Let  us  introduce  into  our  study  the  complex  potential  (z)  of 
a flow  of  an  incompressible  fluid,  equilibrium  at  infinity,  near  a 
certain  foil  2.  Here  represents  the  velocity  at  infinity,  z = x ♦ 
iy  - the  plane  coordinate  of  the  flow  of  incompressible  fluid. 


In  the  neighborhood  of  the  infinitely  distant  point  for  function 
£(z)  - dw/dz  we  have  the  expansion  in  series 


(2.1) 


r® 


C(*)  — 1 + 2V-T.I  * + i*  + • • • • 


where  r*  represents  circulation  over  the  closed  contour  which 
encompasses  foil  l. 


Let  us  use  the  condition 


(2-2) 
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as  a result  of  which  quantity  **  acquires  the  sense  of  coa  plex 
velocity,  quantity  «*  - the  sans*  of  the  velocity  modulus  V of  the 
flow  of  incompreaaibla  fluid.  Osing  this  equality  we  establish  the 
relationship  between  variables  J and  ( and  the  connection  between 
variables  j and  i.  Function  F can  now  be  regarded  as  the  function  of 
the  ccordinate  of  the  plane  of  the  hodograph  i or  as  the  function  of 
the  plane  coordinate  of  the  flew  of  inccu pressible  fluid  z. 

Let  us  determine  function  F(z)  as  follows: 

(2.3)  F(z)  = kVmw(t), 


where  k is  a certain  actual  constant. 


Then,  for  flow  function  'i>  we  get  the  representation 

(2.4)  | = (a), 

while  the  fornula  for  transition  to  the  plane  of  the  gas  flow 
acquires  the  fora  of 

(2.5)  Z=^S*  + ^=[w(*)-M*)]t-- 

with  the  aid  orf  the  last  foraula  we  establish  the  connection 
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between  the  plane  of  the  flow  of  incompressible  fluid  and  the  plane 
of  the  gas  flow.  Here  an  infinitely  remote  point  on  one  plane 
corresponds  to  an  infinitely  remote  point  on  the  other  plane. 
According  to  formulas  (2.2),  (2.5),  and  (1.7)  we  can  plot  in  plane  Z 
a gas  flow  which  is  equilibrium  at  infinity  and  corresponds  to  the 
original  flow  of  incompressible  fluid  - also  equilibrium  at  infinity. 


On  foil  t the  minimal  part  of  the  complex  potential  fMw(s) 
maintains  a constant  value.  If  we  assume  condition  3mFa)u;(i)  = 0 on  t, 
then  on  the  basis  of  (2.4)  and  (2.5)  at  corresponding  points  on 
planes  Z,  determined  according  to  formula 


(2.6) 


Z 


we  will  get  the  equality  't'  * 0.  This  aeans  that  in  plane  Z the  flow 
past  foil  L is  a gas  flew  which  is  uniform  at  infinity,  as  with  the 
original  foil  l . By  the  selection  of  constant  k the  scale  factor  in 
(2.6)  can  be  returned  to  one,  while  foils  t and  L will  be  not  only 
similar  but  also  identical  in  value. 
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Now  let  us  determine  the  connection  between  pressure  coefficient 
Co  at  points  on  the  foil  l i.rd  pressure  coefficient  C ^ at 
corresponding  points  on  L. 

On  the  basis  of  the  third  of  formulas  (1.7)  we  get 


(3.1) 


1 '**  1- 

A2,1*®  (■ 

1 *11 

2n'1  [ 

e2*®  J 

2 »*  | 

Then,  by  using  the  Bernoulli  equation  for  an  incompressible 
fluid,  from  which  it  follows  that 


(3.2) 


C» 1 II 


and  equalities 


(3.3) 


which  emerge  from  formulas  (1.7),  we  arrive  at  the  unknown  connection 


(3.4) 


C — - 

' ✓nwr* 


For  the  connection  v (p)  , determined  by  the  first  and  second  of 


Y..  •* 


DOC 


176  1 


PAGE  9 


equations  (1.7),  the  formula  of  the  linearized  theory  thus  becomes 
precise.  As  a result,  from  formula  (3.4)  we  get,  also  in  the  form  of 
a precise  result,  the  relationship,  known  from  the  linearized  theory, 
between  the  coefficient  of  the  lifting  force  CH  of  a foil  in  a gas 

U 

flow  and  the  coefficient  of  lifting  force  C°  for  the  same  foil  in  a 
flow  of  an  incompressible  fluid: 


(3.5) 


Note  that  for  the  obtained  relationship  v(p)  we  can  give  the 
elementary  proof  of  the  Zhukovskiy  theorem  of  the  lifting  force 
acting  on  a foil.  For  this  we  should  use,  in  place  of  relationships 

(3.5),  (2.6),  and  the  first  and  fourth  from  formulas  (1.7)  , the 

Zhukovskiy  theorem  for  an  incompressible  fluid  and  the  relationship 


(3.6)  r = fcn*(o»+l)r* 

between  circulation  of  velocity  r®  ovmr  the  closed  contour 
encompassed  by  foil  t and  circulation  of  th«  velocity  P over  the 
closed  contour  encompassed  by  foil  L.  Relationship  (3.6)  follows 
directly  from  the  expression  for  the  velocity  pctential  (1.3)  using 
an  expansion  in  series  of  function  £ (z)  . 
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§ 4 


Freedom  in  the  selection  of  constants  n,  B,  and  Ct  can  be  used 
for  approximations  of  the  laws  of  the  adiabatic  flow  of  an  ideal  gas 
or  the  laws  of  other  barotropic  flows  which  interest  us.  Such 
approximations  can  be  plotted  by  different  methods.  Here  we  use  the 
method  discussed  in  [ 1 ],  according  to  which  the  approximation  is  done 
for  the  values  of  parameters  at  a certain  characteristic  point  in  the 
flow.  As  such  a characteristic  point  we  use  an  infinitely  distant 
point,  representing  the  assigned  values  of  gas  parameters  in  it  as 
p , p , v - and  H„. 

roo  uoc 9 oo 9 oo 

The  three  arbitrary  constants  n,  E,  and  Ct  are  obviously 
sufficient  to  satisfy  equalities  v (p^)  = v^,  p (p^)  = pw,  H (p^)  * 

Here  the  equality  K (p^)  = K,,,,  is  satisfied,  where  represents  the 
known  value  of  the  Chaplygin  function  at  the  characteristic  point 
and,  in  addition,  at  an  assigned  point  (p  , p J on  plane  pp  we  have 
the  first  order  of  contact  of  the  assigned  and  approximating  curves 
P (P)  - 


In  Figs.  1*4  dashed  lines  3 represent  approximations  of  the  lavs 
of  the  adiabatic  flow  of  an  ideal  gas  (solid  lines  1)  for  * 0.7. 


The  approximating  dependences  are  determined  by  means  of 
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equation  (1.7)  and  (1.2).  Here 


(4.1) 


n~ 


while  the  foraulas  for  pressure  and  velocity,  after  selection  of 
constants  B and  Ct  , acquire  the  fora  of 


(4.2) 


_p_  , r e»* | v = ««• 

Pea  " 2^/fr-MT  [,2’»  , j ’ v«>  Ooof’®  ’ 


where 


(4.3) 


1 -V1  ~M® 


and  j is  the  adiabatic  exponent.  He  will  assuae  that  quantity  p in 
foraulas  (1.7)  and  (4.1)  is  a dimensionless  quantity  representing  the 
ratio  of  gas  density  to  density  at  the  stagnation  point  of  an 
adiabatic  flow  of  ideal  gas. 


For  coaparison  in  Figs.  1-4  the  curves  plotted  from  the 
approxiaate  Chaplygin  method  (K  = const)  are  shewn  as  thin  solid 
lines  2.  Here  the  approximation  is  done  by  the  same  method  [1]  for 
= 0.7. 


In  both  methods,  for  the  assigned  value  p * p^,  the  Poisson 
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adiabatic  curve  is  approximated  accurate  to  within  small  values  of  (p 
- pj*.  although  within  a broad  range  cf  parameters  change  the 
approximate  Chaplygin  method  has  one  advantage:  Its  laws,  as  we  see 
from  the  curves,  are  in  much  better  agreement  with  the  laws  of  the 
adiabatic  flow  of  ideal  gas. 


Let  us  point  out  certain  conditions  which  limit  the 
applicability  of  this  theory. 

From  the  formula  for  H it  follows  that  <r  - 1/2.  Otherwise,  we 
get  imaginary  values  for  the  M number.  On  the  other  hand,  when  -1  < a 
< 0,  quantity  p,  according  to  the  second  of  formulas  (1.7)  , acquires 
negative  values.  Consequently,  for  physical  considerations,  variable 
• 0 cannot  be  negative:  a } 0.  Since  on  the  basis  of  (2.2)  and  (4.3) 


(5.1) 


a = ln|C|  + 


we  cone  to  the  condition 


(5.2) 


In  I C | > 
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From  this  condition  it  is  specifically  apparent  that  velocity 
should  not  revert  to  zero  at  a single  point  in  the  oriqin&l  flow  of 
the  incompressible  fluid.  This  fact  naturally  imposes  serious 
limitations  on  the  shape  of  the  foil  and  its  position  in  the  flow. 

The  second  limitation  is  related  to  the  formula  for  pressure 

(4.2) .  This  formula  looses  its  physical  sense  at  values  a to  which 
negative  pressure  values  correspond.  From  the  condition  p >y  0 we  get 

(5.3)  Kl<|/ 

Obviously  this  inequality  also  imposes  limitations  on  the  shap>e 
of  the  foil  and  its  position  in  the  flow. 


§ & 


When  0 the  solution  plotted  for  a compressible  flow  moves 

toward  the  solution  for  the  corresponding  problem  of  an 
incompressible  fluid. 

To  be  convinced  of  this,  let  us  transform  the  formulas  for  v and 
Z,  explicitly  isolating  the  dependence  on  parameter  H : 

oo 


(6.1) 


„ _ V ln|C|(l -y/r^M^  + y/rrTilT 

»<*  v°>  yji  - M*, 


(6.2) 


7_,  , i * ~ \Z*  ~~  IP  (z)  — w (») 

_ Z+2  ln|CI(l  - V'  -^)  + V/rri_Wl)  C<*> 


Used  in  the  transformation  process  were  condition  (2.2)  and 
equality  (5.1).  Quantity  k is  determined  frcm  the  condition  of 
equality  to  one  of  the  scale  factors  in  (2.6)  , i.e.,  we  assume 


(6.3) 


Formulas  (6.1)  and  (6.2)  determine  for  an  assigned  the  field 
of  the  velocities  of  the  gas  flow  from  the  known  field  of  flow 
velocities  of  an  incompressible  fluid  (at  corresponding  points  on 
plane  z and  plane  Z the  angles  of  slope  of  the  velocity  vectors 
according  to  condition  (2.2)  are  the  same).  From  (6.1)  and  (6.2)  it 
is  apparent  that  when  11,,-vO,  for  each  fixed  value  z we  get  Z z, 
»/»« V/V^.  Thus,  when  the  field  of  velocities  of  the  gas 

flow  near  the  foil  become  similar  to  the  field  of  velocities  of  the 
original  flow  of  incompressible  fluid. 


It  is  also  simple  to  verify  the  following  transitions  to  the 
limit  shea  p-*  1,  K-+ 1,  r/y. -♦  specifically,  the  last 

transition  to  the  liait  follows  directly  from  relationship 
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(6.4) 


r _ i r* 

v<*  ^1  — MJ,  ' 


which  emerges  after  simple  transformations  from  relationship  (3.6) 
indicated  above. 
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